A Uniform Approximation for the Coherent State Propagator using a Conjugate 
Application of the Bargmann Representation 
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We propose a conjugate application of the Bargmann representation of quantum mechanics. Ap- 
plying the Maslov method to the semiclassical connection formula between the two representations, 
we derive a uniform semiclassical approximation for the coherent state propagator which is finite at 
phase space caustics. 
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Semiclassical approximations have been widely used in 
many areas of physics. They are fundamental to the con- 
ceptual understanding of the quantum-classical connec- 
tion and are also very important in practical situations 
where quantum calculation are difficult, as in systems 
with many of degrees of freedom or with complicated po- 
tential functions. However, semiclassical formulae are of- 
ten not globally valid, i.e., they are not appropriate to de- 
scribe the corresponding quantum function in all regions 
of the space of parameters. The WKB formula for the 
eigenf unctions of a particle in a one-dimensional poten- 
tial well provides a simple example : on the classically 
allowed side of the well the wave-function is oscillatory, 
whereas on the classically forbidden side it is given by a 
single decreasing exponential. At the boundary between 
the two regions, the turning point, the WKB formula 
becomes singular. The formula actually fails in a whole 
neighborhood of the singularity, whose size goes to zero 
as Ti goes to zero. For non-stationary wavefunctions the 
singularities occur at focal points or caustics. After a fo- 
cal point, but sufficiently away from it, the semiclassical 
formulae for wave-functions or propagators still provide 
good approximations, provided the proper Morse phases 
are added @|- 

This general problem of semiclassical expressions, 
which leads to divergences and discontinuities in the 
semiclassical results, can usually be eliminated by prop- 
erly connecting the semiclassical expressions on the dif- 
ferent regions of validity and eliminating spurious con- 
tributions. The most direct way to do that is to solve 
the Schrodinger equation in the vicinity of the singular- 
ity and extend the solution towards the two regions. For 
the WKB problem this amounts to linearize the potential 
about the turning point, leading to the well known solu- 
tion involving the Airy function 0. For non-stationary 
wavefunctions, however, this approach is not usuallypos- 
sible and the Maslov method has to be used 0, 0, 3 . It 
consists basically in changing to a dual representation, 
where the singularity does not exist. For a singularity 
in coordinates, one uses the momentum representation 
and vice- versa. The trick is that, when transforming 
back to the original representation, one should go be- 
yond the quadratic approximation, otherwise the singu- 
larity re-appears. Usually, a stationary phase approxi- 
mation with the exponent expanded up to cubic terms 



is enough, giving rise once again to corrections involving 
Airy functions. 

In this Letter, we study the singularities of the semi- 
classical propagator in the coherent state representa- 
tion. These singularities, called Phase Space Caustics 
(PSCV, have been first identified in [5j and later studied 
in 0,0, 13- All these previous works were concerned with 
ways to identify the singularities and prune the branches 
of spurious contributions arising from them. Here we 
tackle the problem of how to improve the semiclassical 
formula in order to avoid its divergence at the caustics. 
This is a very peculiar situation, since the phase space 
representation provided by the coherent states makes use 
of both coordinate and momentum, leaving no room for 
a natural dual representation. In this paper we define an 
application that works as the canonical conjugate of the 
Bargmann representation 9] and we use it to derive a uni- 
form semiclassical formula for the coherent state propa- 
gator valid in the vicinity of the phase space caustics. For 
the sake of clarity, we restrict ourselves to systems with 
one degree of freedom. Results for multidimensional sys- 
tems, which can be treated along the same lines, and de- 
tailed numerical applications will be published elsewhere. 

The non- normalized coherent state \zq) is defined as 



|zb>=e*» 4T |0>, 
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Here |0) is the ground state of a harmonic oscillator of 
frequency u> = h/mb 2 , is the creation operator and qo, 
Po are the mean values of the position q and momentum 
p operators, respectively. The widths in position b and 
momentum c satisfy be = Ti and zo is complex. The semi- 
classical approximation for the coherent state propagator 
K(z* f ,Zo,T) = (z f \e 



-iHT/h 



zo) is given 



by [Tamoa 
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where M„„ and F depend on (generally complex) classical 
trajectories. These trajectories are best represented in 



2 



terms of the variables u and v, instead of the canonical 
variables q and p, defined by 



1 
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and 
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The sum in Eq. runs over all trajectories governed 
by the complex Hamiltonian H(u,v) = (v\H\u) and sat- 
isfying the boundary conditions u(0) = u' = zq and 
v(T) = v" = z%. Notice that q and p are complex vari- 
ables, while the propagator labels, qo, po for the initial 
state and qf,Pf for the final one, are real. In Eq. (J3J, F 
is given by 

F(v", u', T) = S(v", u', T) + G(v",u\ T) - ~a, (5) 

where S, the complex action of the trajectory, and Q are 
given by 



S(y",u',T) 



in 



(uv - Mil 



H 



, in r ii n i n 
at — — I u v + u v I 



g{v",u',T) = t 



T d 2 H 



dt 



Finally M vv , and its phase a, is an element of the tangent 
matrix defined by 



Su" 
Sv" 



Su' 
Sv' 



(6) 



where Su and Sv are small displacements around the com- 
plex trajectory. We use a single (double) prime to indi- 
cate initial time t = (final time t — T). The elements 
of the tangent matrix can also be written in terms of sec- 
ond derivatives of the action (see Ref. FTa) . Note that 
Eq.JSJ) differs from the formula given in because we 
are using non-normalized coherent states. 

Phase space caustics occur when M TO = 0, causing the 
semiclassical propagator to diverge. Close to these points 
the semiclassical formula provides only a poor approxi- 
mation to the quantum propagator. For a discussion of 
the mechanisms that lead to caustics in systems with one 
degree of freedom see 0, 0, ^| . 

Caustics in the semiclassical propagator in a given rep- 
resentation can usually be circumvented by applying the 
Maslov method. This requires the calculation of the 
semiclassical propagator in the respective conjugate rep- 
resentation, followed by the transformation back to the 
original one, with this last step performed with an ap- 
proximation better than quadratic. For the case of co- 
herent states, although there is no natural dual represen- 
tation to be used with K[v" ,u' ,T), the complex action 
S(v" , u' , T) satisfies the relations 



u{T)= hW 



and 



(7) 



which suggests a Legendre transformation S — * S, by 
the change of variables v" — » (i/Ti)(dS / dv"). The trans- 
formed function S depends on u' and u" , instead of u' 



and v" , 

S{u", u', T) = S(v", u', T) + ihu"v" , (8) 
and satisfies the relations 



v = 



i dS 
hOu 7 * 



and 



i dS 
hdu 1 



(9) 



These properties, on the other hand, suggest the follow- 
ing definition for the dual representation of the semiclas- 
sical propagator: 



K(u",u',T) 



K(v",u',T)e- u " v "dv" (10) 



where the path C will be specified below. 

In the semiclassical limit this integral can be solved by 
the steepest descent method [l5( and an explicit expres- 
sion for K can be obtained. Inserting Eq.© into I|1U|) we 
find the saddle point condition 



± j[s + l nu"v"} = o 



-— (ii) 



where we have considered that Q varies slowly in com- 
parison with S (see Ref. ^3)- Eq.^TJ says that the 
stationary trajectory satisfies u(0) = u' and u(T) = u" , 
i.e., the saddle point value v" of the integration variable 
is equal to v{T) of a trajectory satisfying these bound- 
ary conditions. This imposes that the integration path 
C must coincide with (or be deformable into) a steepest 
descent path through w". Expanding the exponent up to 
second order around this trajectory and performing the 
Gaussian integral we obtain 



K(u",u',T) = , 



|M„ 



S(u" ,u ,T) + j-G(u" ,u ,T)- 



.(12) 



We emphasize that K depends on classical trajectories 
satisfying u' = u(0) and u" = u(T). M uv is given by Eq. 
©, a is its phase, Q(u" , u' , T) is the function Q calculated 
at the new trajectory, and S(u" \u' ,T) is given by Eq. 
©. The PSC affecting K correspond to trajectories for 
which M uv = 0, which generally do not coincide with the 
PSC of K. 

The inverse transformation of Eq. I|1U|) is given by 



K(v",u',T) = 



/,T)e u " v "du" . (13) 



Replacing Eq. (|12fl into i|13|) and doing the integral again 
by the steepest descent approximation up to second order 
gives back the original propagator of Eq. © . 

The pairs of equations ©-1^21) an d (|Tt)jl - ifT3jl look very 
much like the corresponding transformation for the prop- 
agators in coordinates and momenta, K(xf,XQ,T) and 
K(pf,x ,T). In that case the classical trajectory goes 
from xq to Xf in one case and from xq to pf in the 
other. However, although the coherent states define 
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a true quantum mechanical representation [ij, and the 
propagator K{z*^,zq 1 T) corresponds to the matrix el- 
ement (zf\e~ lHT / h \zo) , there is no representation such 
that K(zf, zq, T) also corresponds to a similar matrix el- 
ement. K would involve two kets, \z$) and \zf) instead 
of a ket |zo) and a bra (zf\. 

Therefore, since K is not a matrix element of the evolu- 
tion operator we must formalize its quantum mechanical 
definition so that the previous transformations make pre- 
cise sense. This is done as follows: given a ket |/) and 
its Bargmann representation f(z*) = (z\f) for each 
coherent state ket \w) we define the application 




f(z*)e- z ' w dz*. 

The path 7 must be chosen in such a manner that the 
integral becomes a Laplace transform. This definition is 
suggestive of the need to, so to speak, "cancel the bra 
{z\ and replace it by a ket |to)". At the same time it 
provides just the right Legendre transform we need in the 



semiclassical limit when |/) 
transformation is defined as 



_ p -iHT/h 



/CO = 




f(w)(z\w)dw = 




The inverse 



f(w) e z w dw 



(15) 

with 7' chosen so that the integral is a Mellin transform. 

To illustrate the transformation we apply it to the 
harmonic oscillator. Let |/) be an eigenstate \m) of 
the Hamiltonian operator. Then (z\m) = <j) m {z*) — 



{z*) m /y 



i\ and 



4>m (w) 



V2t 



(z*) m e- z ' w dz* 



(16) 



Writing w = \w\e l6 and z* — re la , the path 7 is defined 
by a = 9 with r varying from to 00. This produces 



4>m (w) 



1 



The inverse transformation is given by 



4>m(z*) 



2ni 



1' 



w 



m+l 



dw. 



(17) 



(18) 



Writing z — \z\e % ^ and choosing 7' so that w = (—a + 
it)e 1 ^ , with a > fixed and t varying from —00 to +00, 
we can solve the integral by the method of residues and 
we find exactly 4> m (z*) = <p m (z*). For the propagator we 
set I/) = e- lHT / h \z ) and obtain 

K(z},z ,T) = (z f \e~ l " T / h \z ) = e z z}e- iuT -i«T/2 (1Q) 

and 

I e -iujT/2 



K{w,z Q ,T) 



z e~ 



(20) 



Equations l|14|l and (|15(l show that the semiclassi- 
cal propagator K given by Ea. (|12J) is the semiclassical 
approximation of a true quantum mechanical function, 
namely Eq.lfT^I with f(z*) = K(z*,z Q ,T). This func- 
tion, although not a matrix element in a mixed repre- 
sentation like K(pj,xo,T) — (j?f\e~ %HT / h \xo), is well de- 
fined provided the integral over 7 converges. The appli- 
cation defined by Eas. l|14|l and (|15l) can be thought of as 
conjugate to the Bargmann representation, and they pro- 
vide the tools to the application of the Maslov method 
to the coherent state propagator. 

The connection between the propagator, Eq.©, and 
its conjugate function, Eq. ((T2*|l . via steepest descent ap- 
proximation with quadratic expansion of the exponent 
works only in the regions where both M uv and M TO are 
non-zero. Close to caustics, where two stationary trajec- 
tories coalesce and M„„ = 0, K is still well defined and 
K can be obtained by doing the inverse transform Q13(l 
using a uniform approximation [l6l | . The basic idea is to 
map the function in the exponent of the integrand into 
an auxiliary cubic function of a new variable X. The new 
function is chosen in such a way that its stationary points 
coincide with those of the original function. Inserting Eq. 
(|12H into Eq. (|13|) we define the new integration variable 
X = X{u") by 

- [S + R] - iu"v" = A - BX + X 3 /3 (21) 

where R = Q + (iln|M u „| — a)h/2 contains the slowly 
varying terms and A and B are functions of u' , v" and 
T. Differentiating both sides with respect to X and dis- 
carding the variation of R yields 



1 dS 
tidu~" 



du" 
~dX 



i[v(T) 



du" 

\~dx 



B + X 2 (22) 



The stationary condition v(T) 



has generally two 



solutions, u" and u"_ , in the vicinity of a caustic. These 
two stationary points, that coalesce at the caustic, are 
mapped into X± — iS 1 / 2 , while the caustic itself cor- 
responds to X — 0. Substituting X = X± in l|21(l and 
solving for A and B we find 



.4 



1 

2h 



R, 



5_ 



B = 



R-) , 

1 2 /3 



R+- R-) 



(23) 



where S± is the action S = S — ihu"v" calculated at the 
stationary trajectories defined by u'±. 

The change of variables from u" to X also produces a 
Jacobian f(X) = du"/dX. Since the X intervals that 
contribute significantly to integral are those close to the 
stationary points, we need to specify the Jacobian only 
in these regions. Writing f(X) = C + G(X - B 1 / 2 ) + 
H(X + B 1 / 2 ) and defining /± = f(X±) wc find 

f(X) = (/+ + /_)/2 + X(f + - f-)/2B^ 2 . (24) 
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FIG. 1: Square modulus of diagonal propagator for fixed z. 
The lines correspond to the exact result (full), bare semiclas- 
sical (dotted) and uniform (dashed). 



be evaluated by the method of steepest descent, with the 
integration path chosen according to the phase of its ar- 
gument z — —B. Therefore the argument of z automat- 
ically indicates whether the two stationary points of the 
exponent contribute to the propagator or if only one of 
them do. Expanding the exponent to second order about 
the contributing points (X + or A_ or both, depending 
on the arg(z)), and doing the resulting Gaussian integral 
recovers the quadratic approximation Eq.©. As an il- 
lustration we consider the Hamiltonian H — (a* a+l/2) 2 . 
Fig.l shows the square modulus of the diagonal propa- 
gator (z\e~ lHT \z) (b = c = H = 1) as a function of T 
for z = 1/(2^/2). The dotted line displays the bare semi- 
classical result, showing a large increase for T > 2.0, due 
to a nearby caustic. The exact result is the full line and 
the uniform approximation is shown by the dashed line, 
which is indeed uniformly good at all times. Detailed 
numerical applications will be published elsewhere. 



Differentiating (|22[1 with respect to X once again and 
calculating at X± we obtain 



f±= ± 



2M ± B 1 ' 2 



1/2 



(25) 



Putting everything together we find the following uni- 
form approximation for the propagator 



K(u',v",T) = 



1 



/2vr 



f{X)S A - BX+x '/*>dX. (26) 



Far from the caustic, where the contribution of each 
stationary trajectory can be computed separately, f(X) 
reduces to f± and the integral can be written in terms 
of an Airy function. In this case the Airy function can 



Eq.J^IJl and the definition of the conjugate application 
and its inverse, Eas. l|14|l and (|15fl . constitute the main 
results of this paper. Although the idea of a conjugate 
application to the Bargmann representation is used here 
just as a tool to derive the above uniform approximation, 
it may be useful in other situations. For instance, a tran- 
sitional approximation, valid only close to the caustics, 
can also be derived. In addition, the Fourier frequencies 
of the transformed propagator are the eigenvalues of the 
Hamiltonian, and it might be simpler to extract those 
eigenvalues from the transformed propagator than from 
the Bargmann propagator. This is certainly the case for 
the Harmonic oscillator, since the time dependence of 
K(w,w,T), see Ea.lfSO)). is trivial. 
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